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ABSTRACT 

In this paper we present a clustering analysis of QSOs over the redshift range z = 0.3 — 
2.9. We use a sample of 10558 QSOs taken from preliminary data release catalogue of 
the 2dF QSO Redshift Survey (2QZ). The two-point redshift-space correlation function 
of QSOs, £q(s), is shown to follow a power law on scales s ~ 1 — 35 hr 1 Mpc. Fitting 
a power law of the form £q(s) = (s/so)~ 7 to the QSO clustering averaged over the 
redshift interval 0.3 < z < 2.9 we find s = 3.99±° 2 3 l h' 1 Mpc and 7 = 1.58±g;Jg for an 
Einstein-de Sitter cosmology. The effect of a significant cosmological constant, Ao, is to 
increase the separation of QSOs, so that with flo = 0.3, Ao = 0.7 the power law extends 
to ~ 60 h- 1 Mpc and the best fit is s = 5.691^ h' 1 Mpc and 7 = 1.56±%$. These 
values, measured at a mean redshift of z — 1.49, are comparable to the clustering of 
local optically selected galaxies. We compare the clustering of 2QZ QSOs to generic 
CDM models with shape parameter r e ff . Standard CDM with T c ff = 0.5 is ruled out 
in both Einstein-de Sitter and cosmological constant dominated cosmologies, where 
r e g ~ 0.2 — 0.4 and r e ff ~ 0.1 — 0.2 respectively are the allowable ranges. 

We measure the evolution of QSO clustering as a function of redshift. For fio = 1 
and Ao = there is no significant evolution in comoving coordinates over the redshift 
range of the 2QZ. QSOs thus have similar clustering properties to local galaxies at 
all redshifts we sample. In the case of £Iq = 0.3 and Ao = 0.7 QSO clustering shows a 
marginal increase at high redshift, sq being a factor of ~ 1.4 higher at z ~ 2.4 than 
at z ~ 0.7. Although the clustering of QSOs is measured on large scales where linear 
theory should apply, the evolution of QSO clustering does not follow the linear theory 
predictions for growth via gravitational instability (rejected at the > 99 per cent 
confidence level). A redshift dependent bias is required to reconcile QSO clustering 
observations with theory. A simple biasing model, in which QSOs have cosmologically 
long lifetimes (or alternatively form in peaks above a constant threshold in the density 
field) is acceptable in an Qq = 1 cosmology, but is only marginally acceptable if 
Qq = 0.3 and Ao = 0.7. Biasing models in which QSOs arc assumed to form over a range 
in redshift, based on the Press-Schechter formalism, are consistent with QSO clustering 
evolution for a minimum halo mass of ~ 10 12 M© and ~ 10 13 M Q in an Einstein-de 
Sitter and cosmological constant dominated universe, respectively. However, until an 
accurate, physically motivated, model of QSO formation and evolution is developed, 
we should be cautious in interpreting the fits to these biasing models. 

Key words: galaxies: clustering - quasars: general - cosmology: observations - large- 
scale structure of Universe. 



1 INTRODUCTION 
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The 2dF QSO Redshift Survey (2QZ) aims to compile a 
homogeneous catalogue of ~ 25000 QSOs using the Anglo- 
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Australian Telescope (AAT) 2-degree Field facility (2dF; 
Taylor, Cannon & Watson 1997). This catalogue will con- 
stitute a factor of J>50 increase in numbers to a equivalent 
flux limit over previous data sets (e.g. Boyle et al. 1990). 
The main science goal of the 2QZ is to use QSOs to probe 
the large-scale structure of the Universe over a range of 
scales from 1 to 1000 h~ x Mpc, and in the redshift inter- 
val, 0.3<z<2.9. 

Clustering of QSOs at small to intermediate scales 
(1 — 50 /i" 1 Mpc) supplies a wealth of information on large- 
scale structure. QSOs still give us the only method of di- 
rectly determining the 3-dimensional clustering of high red- 
shift objects within a large enough volume for it to be truly 
representative. When complete, the 2QZ will sample a vol- 
ume of 1.5 x 10 9 h _3 Mpc 3 (for f2o = 1), an order of mag- 
nitude larger than current galaxy redshift surveys (e.g. the 
2dF Galaxy Redshift Survey; Colless 1999). This large vol- 
ume also allows us to probe the scales where linear evolution 
occurs, simplifying comparisons with theory. 

The shape and amplitude of the two-point auto- 
correlation function, £(r), is determined by two factors. The 
first is the distribution of matter fluctuations in the Uni- 
verse. This depends on fundamental physics, such as the 
growth of structure via gravitational instability and the ini- 
tial spectrum of fluctuations. The second factor concerns the 
complex and generally non-linear physics which occurs dur- 
ing galaxy and QSO formation. The difference between the 
matter and galaxy or QSO distributions is commonly called 
bias, b(r,z), such that 



(1) 



where £cj(r, z) and £ p (r,z) are the two-point correlation 
functions of QSOs and the density field respectively. Both 
are functions of scale, r, and redshift, z. Often, a linear bias 
is assumed, which has no scale dependence, and it appears 
likely that for any local process of galaxy formation b should 
tend to a constant value on scales where the density per- 
turbations are linear (e.g. Mann, Peacock & Heavens 1998; 
Peacock 1997). We will assume a linear bias throughout this 
paper. 

The first attempt to measure the clustering of QSOs 
was made by Osmer (1981). Shaver (1984) was the first to 
detect QSO clustering on small scales, although in an in- 
homogeneous sample. Shanks et al. (1987) made the first 
detection of clustering at <^10 ft -1 Mpc in a complete and 
uniform ly selected sample ; part of the Durham/AAT UVX 
survey (Boyle et al. 199C). A number of authors have used 
this and other QSO samples to measure clustering. They all 
reach generally the same conclusions that clustering is de- 
tected at the ~ 3 — 4<r level and is approximately consistent 
with a clustering scale length ro ~ 6 h^ 1 Mpc, simi lar to lo- 
cal gal axy clustering, at a mean redshift of z ~ 1.4 ( lovino & 



Shaver |l988fe |Andreani fc Cristiani 1992|; |Mo fc Fang 1993f 



Shanks & Boyle 1994; Croom & Shanks 1996). There has 



been significant disagreement over the redshift evolution of 
QSO clustering including claims for a decrease in the QSO 
correlation length (ro) with redshift (lovino & Shaver 1988), 



an increase in ro with redshift (La Franca, Andreani & Cris- 
tiani 1998; hereforth LAC98) and no change with redshift 
(Croom & Shanks 1996; hereforth CS96). 

The measurement of galaxy clustering at high redshift 
has also taken dramatic steps forward in recent years. A 



number of surveys have made measurements of the cluster- 
ing strength of galaxies up to z ~ 1. These samples typically 
contain a few hundred to a thousand galaxies over relatively 
small areas. The Canada-France Redshift Survey (CFRS) 
shows a signific ant decrease in cluster ing amplitude on scales 
< 2 h^ 1 Mpc ( Le Fevre et al. 1996[). However, larger sam- 
ples, such as the CNOC-2 survey ( parlberg et al. 199£ ) show 
much slower evolution, with a gradual decrease of clustering 
with redshift: r (z) oc (l + z)~°- 3±0 ' 2 . Deep wide-field (~ few 
degrees) imaging surveys used to measure the angular cor- 
relation function of galaxies also sug gest higher clusterin g 
amplitudes than found in the CFRS (Postman et al. 1998). 
The differences found between these samples is partly due 
to the different selection methods (e.g. magnitude limits; 
photometric bands) used, ft is likely that different galaxy 
types cluster differently, e.g. optically vs. infrared selected 
galaxies (Peacock 1997). Clustering is also likely to be a 
function of galaxy luminosity. It is possible that this is the 
case for QSOs, although we leave the discussion of lumi- 
nosity dependent clustering of QSOs to a future paper. We 
note, however, that due to the stong luminosity evolution of 
QSOs (e.g. Boyle et al. 2000) an apparent magnitude lim- 
ited survey of QSOs samples approximately the same part 
of the luminosity function at all redshifts up to z ~ 2. A 
second affect responsible for the difference in galaxy clus- 
tering results is cosmic variance due to the small volumes 
and scales sampled, in particular by the CFRS. A key ele- 
ment of the 2QZ is that it is large enough to minimize any 
effects of cosmic variance on scales smaller than a few hun- 
dred Mpc. Although studies of galaxy clustering have 
been typically limited to z<l, Steidel et al. (1998) have used 
galaxies detected by their Lyman-break to derive the clus- 
tering properties of galaxies at z ~ 3. These observations 
show that the clustering of L ~ L* galaxies at z ~ 3 is 
also similar to local galaxies on scales <10 Mpc, with 
ro ~ 4 — 6 h~ x Mpc de pending on the assumed cosmology 



( |Adelberger et al. 1998|) . 

In this paper we look at QSO clustering in the 2QZ 
on scales from ~ 1 to 100 h^ 1 Mpc. We do not attempt 
to study larger scales because of the current non-uniformity 
of the data set. This will be reserved for future work, on 
completion of the survey. In Section ^ we describe the 2QZ 
data used and our methods of analysis. In Sections ^ and 
^ we present our clustering results and compare them to 
physical models. Our conclusions are presented in Section || 



2 DATA AND TECHNIQUES 

2.1 The 2dF QSO Redshift Survey 

For the analysis in this paper we have used the first public re- 
lease catalogue of the 2QZ containing 10681 QSOs (the 1 0k 
catalogue). This 10k catalogue contains the most spectro- 
scopically complete fields observed prior to November 2000 
and will be released to the astronomical community in the 
first half of 2001. The sample contains 10558 QSOs in the 
redshift range 0.3 < z < 2.9 which will be included in our 
analysis below. 

The identification of QSO candidates for the 2QZ was 
based on broad band ubjr colours from Automatic Plate 
Measuring (APM) facility measurements of UK Schmidt 
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Figure 1. The current fractional coverage in the NGC (top) and SGC (bottom) strips of the 2QZ. Each circular region corresponds to 
one 2dF pointing. The coverage is the fraction of QSO candidates observed in each region. The small rectangular holes correspond to 
regions containing bright stars and plate defects. 



Telescope (UKST) photographic plates. The survey com- 
prises 30 UKST fields, arranged in two 75° x 5° declination 
strips centred in the South Galactic Cap (SGC) at S — —30° 
and the North Galactic Cap (NGC) at S = 0° with RA 
ranges a = 21 h 40 to 3 h 15 and a = 9 h 50 to 14 h 50 respec- 
tively .F^ach_TJKSTJidd con^ CCD calibra- 



tion (Boyle et al. 1995; Croom et al. 1999). The completed 



survey will cover approximately 740 deg (some areas having 
been removed due to bright stars, plate defects etc). Further 
details of the photometric catalogue can be found in Croom 
(1997), Smith (1998) and Smith et al. (2001). 

Spectroscopic observations have been carried out using 
the 2dF instrument at the AAT in conj unction with the 
2dF Galaxy Redshift Survey ( |Colless 1999b , as the 2QZ and 
galaxy survey areas cover the same regions of sky. Typically, 
225 fibres are devoted to galaxies, 125 to QSO candidates 
and 25-30 to sky in each 2dF observation. Spectroscopic data 



are red uced using the 2dF pipeline reduction system (Bai- 



ley & Olazebrook 199!:). The identification of QSO spectra 
and redshift estimation was carried out using the AUTOZ 
code written specifically for this project (Miller et al. 2001 
in preparation). This program compares template spectra of 
QSOs, stars and galaxies to the observed spectra. Identifica- 
tions are then confirmed by eye for all spectra. Spectroscopic 
completeness is typically > 80 per cent when observations 
are made in reasonable or good conditions. In the analysis 
below we use all objects which have been classified as class 
1 QSOs (class 1 being the highest quality identification; ob- 
jects classified as class 2 IDs or "QSO?" were not included) 
and which were observed in fields within the 10k catalogue 
(which is limited to > 85 per cent spectroscopic complete- 
ness) . 

2.2 Correlation function estimates 

As the QSO correlation function, £q probes high redshifts 
and large scales, the measured values are highly depen- 
dent on the assumed cosmology. We employ the method of 
Osmer (1981), which uses the coordinate transform in the 
Robertson- Walker metric (Weinberg 1972) to determine the 



comoving separation of pairs of QSOs. We choose to calcu- 
late £q for two representative cosmological models; f2o,Ao = 
(1,0) and (0.3,0.7), where fio and Ao represent the con- 
ventional mass and vacuum energy (cosmological constant) 



density contribution respectively, to the total energy density 
of the Universe. We will call these cosmological models EdS 
(Einstein-de Sitter) and A respectively. 

We have used the minimum variance estimator sug- 
gested by Landy & Szalay (1993) to calculate £(s), where 
s is the redshift-space separation of two QSOs. This estima- 
tor is 



QQ(s) - 2QR(s) + RR(s) 
RRjs) 



(2) 



where QQ, QR and RR are the number of QSO-QSO, QSO- 
random and random-random pairs counted at separation s± 
As. We bin our pairs such that log(As) = 0.1 or 0.2. The 
density of random points used was 50 times the density of 
QSOs. 

The area of the survey is covered by a mosaic of 2dF 
pointings. These pointings overlap in order to obtain com- 
plete coverage in all areas, including regions of high galaxy 
and QSO density. As the survey is not yet complete this 
means that certain areas within 2dF fields will not have 
had all candidates observed, and therefore the observational 
completeness of the sample varies strongly with angular po- 
sition on the sky. This variation in observational complete- 
ness can clearly be seen in Fig. |l| Where a large number of 
2dF pointing overlap the coverage is ~ 100 per cent, while 
in overlap regions which have yet to be observed a second 
or third time the completeness is significantly lower. Par- 
ticular care has been taken to construct the random point 
distribution so as to take into account this angular selection 
function. In each region defined by the intersection of 2dF 
fields we have counted the number of QSO candidates ob- 
served and compared this to the total number to calculate 
the fractional observational completeness. We then weight 
the probability of a random being placed in that region by 
this fractional completeness. This corrects for the angular 
incompleteness due to overlapping 2dF fields. 

Our candidate density is not completely uniform over 
the length of the strips, due to an increase in stellar con- 
tamination in areas closer to the galactic plane. Secondly, 
small residual calibration errors in the relative magnitude 
zero points of the UKST plates could add spurious struc- 
ture on large scales. Any possible offsets are being corrected 
by calibration from further CCD photometry, however in 
this paper we will correct for this effect by normalizing the 
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Figure 2. The redshift distribution of QSOs in the 10k cata- 
logue for the NGC (dotted histogram) and SGC (solid histogram) 
strips. Also shown is the fit used to generate random redshift dis- 
tributions (smooth curves), normalized to the number of QSOs 
in each slice. 



number of random points to the number of QSOs with spec- 
troscopically determined redshifts in each UKST field. This 
correction will clearly remove power on large scales, which 
is why we do not discuss structure on scales larger than 
~ 100 hT 1 Mpc in this paper. After constructing the angular 
mask, we then assign the random points a random redshift, 
taken from a spline fit to the binned (Az = 0.2) redshift 
distribution of the full 2QZ sample (changing Az by a fac- 
tor of 2 makes no observable difference to £q). The redshift 
distributions for the NGC and SGC data sets are shown in 
Fig [| The fit used to generate the redshift distribution for 
the random points is also shown (the smooth curves in Fig. 
^|), in each case normalized to the number of QSOs in each 
2QZ strip. The above process for correcting observational 
completeness and calculating £q we call method 1. 

We have tested the effectiveness of this process by mak- 
ing comparisons to correlation functions derived using two 
other methods. The first (method 2) is to calculate the cor- 
relation function from regions of the survey that have no 
overlapping fields still to be observed, that is, they have 
~ 100 per cent observational coverage. The number of QSOs 
in these regions is significantly less that in the total sample, 
reducing the signal-to-noise in £q. The second comparison 
method (method 3) is to allocate each random point an (a, 5) 
taken from the QSO catalogue, so that the random distribu- 
tion has exactly the same angular distribution as the QSOs. 
The redshifts of the random points are then allocated using 
the spline fit discussed above. 

Two other observational biases could, in principle, af- 
fect our measurements of £q. The first is due to the fact 
that the 2dF instrument cannot position two fibres closer 
than ~ 30". We are therefore currently biased against small 
angular separation QSO pairs (this problem is being reme- 
died by independent follow-up of close QSO pairs). We have 
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Figure 3. The angular correlation function for all currently ob- 
served candidates in the 2QZ, split into the NGC (filled circles) 
and SGC (open circles). A deficit of pairs on small scales due to 
2dF positioning constraints can be seen. The dotted line denotes 
the weight function used to correct for the lack of close pairs. The 
uniformity at large scales, > 0.1°, demonstrates the effectiveness 
of our correction for the non-uniform field coverage. 



measured the angular correlation function of observed can- 
didates, which shows this bias (see Fig. ^). Measuring the 
extent of the anti-correlation in Fig. ^ allows us the correct 
for the close pairs bias. The dotted line, which traces the 
anti-correlation is uj (6) = 4.0 x 10~ 5 6i- 2 , and this can be 
used to construct a function 



W cp (6) = 



1 



1 - 4.0 x 10- 5 6>- 2 : 



(3) 



which is the weight function for close pairs separated by 9 
degrees. In practice this correction makes no difference to 
the measured correlation function as almost all of the QSO 
pairs with small angular separations have widely differing 
redshifts, and the weighting of a small number of pairs has 
a negligible effect on large scales. 

Extinction by galactic dust will also imprint a signal on 
the angular distribution of the QSOs. Primarily this changes 
the effective magnitude limit in bj by Abj = 4.035 xE(B-V) 
where we use the dust reddening E(B — V) as a function of 
position calculated by Schlegel, Finkbeiner & Davis (1998). 
We then weight the random distribution according to the 
reduction in number density caused by the extinction such 
that 



W cxt (a,8) = W~ 0A ^ (o 



(4) 



where (3 is the slope of the QSO number counts at the mag- 
nitude limit of the survey. At bj — 20.85, the magnitude 
limit of the 2QZ, the QSO number counts are flat, with 
/3 ~ 0.3. Again we find that applying this correction makes 
no significant difference to the measured £q. 

It can be useful to present clustering results in a 
non-parametric form, specified by the clustering amplitude 
within a given comoving radius, rather than as a scale length 
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Figure 4. The distribution of 2QZ QSOs in the 10k catalogue. Note that only objects at 0.3 < z < 2.9 are used in our analysis. The 
SGC strip is on the left, the NGC on the right. The rectangular regions show the distribution projected onto the sky in each strip. 



which depends on a power law fit to £q. This is generally 
represented by the integrated correlation function, £, within 
a given radius in redshift-space, s max , 



£,{x)x Ax. 



(•») 



Authors tend to choose a variety of values fo r s 

s m ax = 10 h" 1 Mpc ( |Shanks fc Boyle 1994 |Croom & 
Shanks 



1996) or s n 



15 h Mpc (La Franca et al 



1996). The choice is a compromise, selecting the scale for 
which a significant signal is seen. It is easiest to relate these 
measurements to theory for large scales, where linear evolu- 
tion occurs. Below we will quote clustering amplitudes with 
Smax = 20 ft -1 Mpc as this is a scale at which evolution 
should be linear to better than a few per cent. We note that 
choosing a large radius also reduces the effects of small scale 
peculiar velocities and redshift measurement errors, which 
may well be a function of redshift. 

We calculate the errors on £q using the Poisson estimate 

of 



A£(s) = 



VQQ{s) 



(6) 



At small scales, <J50 ft -1 Mpc, this estimate is accurate be- 
cause each QSO pair is independent (i.e. the QSOs are not 
generally part of another pair at scales smaller than this). On 
larger scales the QSOs pairs become more correlated and we 
use the approximation that A£(s) = [1 + £(s)]/ WNq, where 
Nq is the total number of QSOs used in the analysis (Shanks 



& Boyle 1994; CS96). In this paper, we will generally be con- 
cerned with analysis on small scales (< 50 ft -1 Mpc), where 
the Poisson error estimates are applicable. As a confirmation 
of our Poisson error estimates we have also derived field-to- 
field errors, by splitting the NGC and SGC strips into two, 
and determining the scatter between the resulting four in- 
dependent regions. The errors determined in this fashion 
are approximately equal to or less than the Poisson errors. 
We also test bootstrap errors which are found to be ~ V3 
times greater than Poisson on all scales of interest, in agree- 
ment with expected theory (Mo, Jing & Borner 1992) and 
previous measurements (e.g. Boyle & Shanks 1994; CS96). 
On small scales, <C2 hT 1 Mpc, the number of QSO-QSO 
pairs can be <10. In this case simple root-n errors (Eq. ^|) 
do not give the correct upper and lower confidence limits 
for a Poisson distribution. We use the formulae of Gehrels 
(1986) to estimate the Poisson confidence intervals for one- 
sided 84% upper and lower bounds (corresponding to lcr for 
Gaussian statistics) . These errors are applied to our data for 
QQ(s) < 20. By this point root-n errors adequately describe 
the Poisson distribution. 



2.3 Fitting models to £(s) 

Below we make comparisons of the data to a number of 
models, both simple functional forms (power laws) and more 
complex, physically motivated, models (e.g. CDM). We use 
the maximum likelihood method to determine the best fit 
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Figure 5. The two-point correlation function for 2QZ QSOs in the redshift interval 0.3 < z < 2.9 for an EdS cosmology, using the 
different estimators discussed in the text, a) method 1, b) method 2, c) method 3 and d) all three methods. In plots a), b) and c) we 
show £q from the combined NGC and SGC strips with the best fit power law at s < 35 h^ 1 Mpc in each case. In d) we compare the 
three methods, the two points for each method being separate estimates in the NGC (filled symbols) and SGC (open symbols) strips. 
The lines are identical to those in a), b) and c) denoting the best fit power law in each method. 



parameters. The likelihood estimator is based on the Poisson 
probability distribution function, so that 

i-n^ ^ 

is the likelihood, where v is the observed number of QSO- 
QSO pairs, fi is the expectation value for a given model and 
N is the number of bins fitted. We fit the data with bins 
Alog(r) = 0.1, although we note that varying the bin size 
by a factor of two makes no noticeable difference to the resul- 
tant fit. In practice we minimize the function S = — 21n(L), 
and determine the errors from the distribution of AS, where 



AS is assumed to be distributed as x 2 . This procedure does 
not give us an absolute measurement of the goodness-of-fit 
for a particular model. We therefore also derive a value of 
X 2 for each model fit in order to confirm that it is a reason- 
able description of the data. In particular this is appropri- 
ate when fitting on moderate to large scales (J>5 hT 1 Mpc), 
where the pair counts are large enough that the Poisson er- 
rors are well described by Gaussian statistics. 
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3 THE CORRELATION FUNCTION OF 2QZ 
QSOS 

Here we present the results of our clustering analysis on an 
initial sample of 2QZ QSOs. This sample contains 10558 
QSOs taken from the 2QZ 10k catalogue. Fig. ^ shows the 
distribution of QSOs projected onto a plane of constant dec- 
lination. We note that the current distribution is highly non- 
uniform, as the survey is only partially complete. 



3.1 The redshift averaged QSO correlation 
function 

We first measure the QSO two-point correlation function av- 
eraged over the entire redshift interval 0.3 < z < 2.9. For 
an EdS cosmology we estimate £q using the three differ- 
ent processes discussed in Section 2.2: 1) full accounting for 



non-uniform coverage, 2) taking only completely observed 
regions, and 3) using the QSO (a,S) for the random point 
positions. These are presented in Fig. ^a, b and c respec- 
tively. The results demonstrate that QSO clustering follows 
a power law on small to intermediate scales. There is some 
evidence of a break in the power law at ~ 35 h~ l Mpc. We 
fit a power law of the conventional form, 



as) 



(8) 



The best fits using the maximum likelihood technique are 
(50,7) = (3.99l°;g,1.58±°;i»), (4.59t^, 1.64±g;g) and 
(3.87±g;l|, 1.63±g;^) for methods 1, 2 and 3 respectively, 
where so is in units of h' 1 Mpc. We fit the power law on 
scales 0.7 — 35 h" 1 Mpc. The minimum scale is set by the 
smallest scale at which we find QSO pairs and the maxi- 
mum scale is set by scale of the observed break in £(s). A 
comparison of all three methods is shown in Fig. ^d. Here 
we also plot separately the clustering of the NGC and SGC 
strips. 

First, we note that the signals from the NGC and SGC 
strips are consistent. The NGC has no pairs at very small 
scales (< 1.5 ft -1 Mpc), however the SGC strip only contains 
3 pairs at these scales, and fewer are expected in the NGC 
due to the smaller number of objects in this strip (4005 in 
the NGC vs. 6553 in the SGC). Second, there appears to be 
no significant difference between our different estimations of 
£q. Method 2 shows a slightly higher signal while method 
3 is marginally lower than the other two methods. We esti- 
mate how much of the difference between methods 1 and 3 
could be due to the removal of real signal by taking the mea- 
sured correlation function from method 1 and integrating it 
over our redshift range, weighted by the QSO redshift dis- 
tribution. This then gives us an angular correlation function 
with which we weight the random distribution when deriving 
the 3-D correlation so suppressing the angular component. 
The results of this analysis are shown in Fig. ^| On small 
scales there is very little effect on the 3-D clustering, how- 
ever on scales >20 ft -1 Mpc the clustering signal becomes 
suppressed by larger amounts (dotted line). The correlation 
functions measured from the data using the two methods 
are also plotted in Fig. []. The difference in the measured 
values at J>20 Mpc is similar to that predicted by the 
model, suggesting that some large-scale power is removed 
by method 3. We therefore choose to use method 1 through- 
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Figure 6. The estimated effect of using actual QSO (»,<5) for 
the random distribution when estimating the correlation function 
(method 3; dotted line). The solid line is the input model power 
law (identical to the power law fit to method 1). The filled and 
open circles are the estimates of §q from the data using methods 
1 and 3 respectively. 
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Figure 7. A comparison between the 2QZ §Q for different cos- 
mologies: EdS (filled circles) and A (open circles). The solid and 
dotted lines are the best fit power laws in each case. 



out the remainder of our analysis (method 2 contains half as 
many QSOs as method 1, only 5348 and they are generally 
distributed in many small overlap regions; the dark shaded 
regions in fig. |l|). Any residual systematic errors caused by 
the variable observational completeness are not significant 
enough to affect any of the conclusions of this paper. 
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Figure 8. The two-point correlation function for 2QZ QSOs in the redshift interval 0.3 < z < 2.9 in an EdS cosmology, compared to 
the clustering of local galaxies from the Durham/UKST Survey (Ratcliffe et al. 1998, triangles) and the Las Campanas Survey (Tucker 
et al. 1997, squares), plotted in a) log-log space to highlight smaller scales and b) log-linear space to highlight larger scales. In a) the 
2QZ data is plotted in log(As) = 0.2 bins, while in b) log(As) = 0.1. The solid line is the best fit power law to the 2QZ data for the 
EdS cosmology, while the dotted line is the best fit for the A cosmology. The 2QZ data points for the A model are omitted for clarity. 



Table 1. 2QZ clustering results for various cosmologies and redshift intervals. The so an d 7 are best fit values. The results for the 2 
parameter fit are allowing both so and 7 to vary freely. For the 1 parameter fit we constrain 7 to be the best fit value for each cosmology 
over the full redshift interval (0.30 < z < 2.90) and allow only sq to vary. The reduced x 2 f° r each fit are also listed. 



(fio.Ao) 


redshift 


range 


z 




2 parameter fit 

so 7 


x 2 


1 paramet 
so 


sr fit 

x 2 


£(20) 


(1.0,0.0) 
(0.3,0.7) 


0.30 < 
0.30 < 


z 
z 


< 2.90 

< 2.90 


1.49 
1.49 


10558 
10558 


3 qq+ - 28 
5 69+ - 42 


1 50+O.IO 
i -°°-0.09 
1 5R+O.IO 
1 ' OD -0.09 


1.42 
1.33 






0.197 ±0.026 
0.416 ±0.048 


(1.0,0.0) 
(1.0,0.0) 
(1.0,0.0) 
(1.0,0.0) 
(1.0,0.0) 


0.30 < 
0.95 < 
1.35 < 
1.70 < 
2.10 < 


z 
z 
z 
z 
z 


< 0.95 

< 1.35 

< 1.70 

< 2.10 

< 2.90 


0.69 
1.16 
1.53 
1.89 
2.36 


2299 
2116 
2177 
2186 
1780 


3 o 4 +0.56 

9 79 +0.94 
Z - '^-1.18 
, 4q +0.6l 

4 Ql+0.55 

-0.61 
4 43+0.77 
^•*°-0.94 


1 70+ ' 27 
!• ,u _0.36 
1 9C+0.27 
i '^°-0.25 

1 63+ 022 
u -0.22 
1 00+O.21 
x.oo_q 20 

1 84+0.30 


1.14 
1.37 
1.51 
0.64 
1.14 


q fi5 +0.56 
°- DO -0.56 
q 55+O.6I 
• 3 - OO -0.64 
3 41+0.56 
°'^ 1 -0.56 

q O5+0.56 

3 QR+0-80 
°- ao -0.83 


1.18 
1.27 
1.53 
0.87 
1.35 


0.163 ±0.054 
0.211 ±0.057 
0.192 ±0.052 
0.140 ±0.055 
0.099 ± 0.078 


(0.3,0.7) 
(0.3,0.7) 
(0.3,0.7) 
(0.3,0.7) 
(0.3,0.7) 


0.30 < 
0.95 < 
1.35 < 
1.70 < 
2.10 < 


z 
z 
z 
z 
z 


< 0.95 

< 1.35 

< 1.70 

< 2.10 

< 2.90 


0.69 
1.16 
1.53 
1.89 
2.36 


2299 
2116 
2177 
2186 
1780 


r 9Q+0.72 
°- ZO -0.89 

t- 9 o+o.92 

J - zo -1.08 

6 24+ - 86 
D - Z -0.99 
R q3+l-32 
D - M,3 -1.64 


1 7 9+0.23 
1 -' z -0.22 
-, qo+0.27 
1-3S_ 24 

1 55+0.21 
i.oo_ jo 


1.21 
0.69 
1.92 
0.83 
1.26 


4 90+ ' 71 

4 65+ - 89 
^•DO-0.91 

5 24+ - 82 
°- z ^-0.81 
5 54+0-84 


1.28 
0.63 
1.92 
1.09 
1.30 


0.269 ±0.085 
0.371 ±0.102 
0.468 ±0.103 
0.394 ±0.110 
0.615 ±0.178 



We also calculate £q for the A cosmology, with Qq = 0.3 
and Ao = 0.7. This is compared to the method 1 estimate 
for the EdS case in Fig The effect of introducing a sig- 
nificant cosmological constant term is to increase the rela- 
tive separation of QSOs, and hence increase the clustering 
scale length. The break in the power law is now seen at 
~ 60 h' 1 Mpc, we therefore make our power law fits out 
to this scale. The best fit power law for the A cosmology is 



s , 7) = (5.69to;5o, 1.56±S;Si). AH results are listed in Table 



3.2 QSO clustering compared to local galaxies 

In Fig. ^ we compare our QSO clustering results at z = 1.49 
to galaxy clustering at low redshift (z ~ 0.05). In particular 
t he Las Campanas ([T ucker et al. 1997) and Durham/UKST 



(Ratcliffe et al. 19981) galaxy surveys (open squares and tri- 
angles respectively). We see that there is good general agree- 
ment between the galaxy and QSO clustering, although the 
samples have differing redshift ranges. The EdS £q is slightly 
lower on average than £ ga i, while £q in the A cosmology is 
closer in amplitude to the galaxies. Both QSOs and galaxies 
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Figure 10. A comparison between the 2QZ §q and and CDM type models in the a) EdS and b) A cosmologies. The dotted lines are 
non-linear CDM correlation functions with various values of F e g scaled by a linear bias to give the best fit to the 2QZ data. r c ff = 0.1, 
0.2, 0.3, 0.4 and 0.5 from top to bottom at large scales). 
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Figure 9. A comparison between the 2QZ £q (open circles) and 
£ p (dotted line) from the Hubble Volume simulations in the A 
cosmology. The dashed line is 2.1 2 X £ p . The arrow marks the 
resolution limit of the simulation. Also shown is the best fit power 
law (solid line). 



show a break in £(r) at ~ 40 h 1 Mpc. We note that the 
errors on £q are smaller than those on £ ga i at > 20 h^ 1 Mpc. 



3.3 QSO clustering compared to CDM 

In order to compare directly to theory, and include all non- 
linear effects and redshift-space distortions, we have used 
the Hubble Volume simulations of the Virgo Consortium 



( |Colberg et al. 199oj ). We have produced mock 2QZ QSO 
catalogues with the same survey geometry and explicitly in- 
cluded evolution of the density field by outputting the sim- 
ulation at different times along the light cone. A detailed 
discussion of the simulations, including a number of bias- 
ing models will be given in Hoyle et al. (2001 in prepara- 
tion). Here we simply compare the dark matter correlation 
function averaged over the light cone to the 2QZ data. In 
particular we compare the redshift-space mass correlation 
function of a ACDM model to £q. This model has Qo = 0.3, 
Obaryon = 0.04, Ao = 0.7, erg = 0.9, h = 0.7 and an effective 
shape parameter, tak ing into account the baryon compo- 
nent ( Bugiyama 1995 ), of r e ff = 0.17. The model and data 
are shown in Fig. |^. The amplitude of £q is ~ 4 times larger 
than the ACDM mass correlation function, £ p . When scaled 
by this factor, the model and data appear to be well matched 
with a best fit bias value of 2.1. 

Increasing the value of r e fi will move the models away 
from the data by steepening £ p at large scales. We do not 
have a large suite of simulations with which to compare 
the effect of changing r c ff and cosmology. However, on the 
scales which we are fitting, linear theory is a reasonable ap- 
proximation. Therefore the effect of redshift space distor- 
tions will be simply to sc ale £ by (1 + 2/3/3 + /3 2 /5) where 



f3 ~ fi(]' 6 /6 (Kaiser 1987). We can then simply absorb this 



factor into an effective linear bias factor. We then fit model 
real space non-linear correlation functions at z — 1.49 to the 
data (again at scales 5 to 100 ft -1 Mpc) using the ansatz of 
Peacock & Dodds (1996) to determine the non- linear cor- 
rection to the model £. The deviation from non-linearity is 
small (typically <J5%) on the scales of interest. We do not 
take into account small-scale non-linear velocity dispersions 
in our model, however these should be small at the scales 
and redshifts considered. We also do not consider the effects 
of redshift measurement errors on £(s), these again should 
only be a factor on small, <5 h~ x Mpc, scales. We use five 
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different models with r eff = 0.1, 0.2, 0.3, 0.4, 0.5 and fit for 
the effective linear bias value. The results of this procedure 
are shown in Fig. In the EdS cosmology, models with 
r e ff = 0.2, 0.3 and 0.4 are acceptable at the 10% level, while 
r c ff = 0.1 and 0.5 are ruled out at greater than 90% con- 
fidence. The main reason that a broad range of models are 
acceptable is the relatively low point at ~ 20 h -1 Mpc in 
£q. In the A cosmology the r e ff =0.1 and 0.2 models are 
the only ones to agree with the data, the others being ruled 
out at greater than 99.9% confidence. Thus the QSO corre- 
lation function detects excess large-scale power over what is 
expected in the V e g = 0.5 standard CDM model. Confirm- 



ing the results from the APM galaxy survey ( Maddox et al 
199(j)~| 

The required r e fi is larger in the A cosmology, as struc- 
ture is moved to larger scales. This suggests a test with the 
full 2QZ which will be devoid of observational incomplete- 
ness as well as having increased statistical accuracy. The 
break of the correlation function in a CDM type cosmology 
can be used as a standard rod to determine cosmological 
parameters, in particular Ao, if it is at linear scales. For ex- 
ample, if at low redshift the shape is well defined, then if 
the break is in the linear theory regime it should remain at 
the same scale at high redshift. Measuring the break at a 
different scale at high redshift would imply the wrong cos- 
mological parameters were being used in the determination 
of the high redshift correlation function. This is similar to 
the geo metric tests discussed by several authors ( Alcock & 
Paczyn jjki 197S| ; [Phillipps 1994] ; jBallinger et al. 1996| ) but has 
the advantage of not being affected by redshift-space distor- 
tions if clustering can be measured on a sufficiently large 
scale. This is because linear redshift space distortions only 
affect the amplitude and not the shape of £. Shanks & Boyle 
(1994) proposed a similar method, using linear features in 
the correlation function on <;100 h~ Mpc scales. 



4 EVOLUTION OF QSO CLUSTERING 

4.1 Measurements of QSO clustering evolution 

In the previous section we calculated £q averaged over a 
large redshift interval. We now split the 2QZ QSO sample 
up into five redshift intervals containing approximately equal 
numbers of QSOs. The exact limits and numbers of QSOs 
are given in Table [lj. The measured £q are shown in Fig. 
pd] for the EdS cosmology. QSO clustering appears to vary 
little over the entire redshift range we consider. The data 
points are consistent with the redshift averaged £q (dotted 
line in Fig. [ll]). For each redshift interval we fit a power 
law, the results of which are shown by the solid lines in Fig. 
pd] and in Table [lj As for the redshift averaged analysis we 
fit the power law within 35 h^ 1 Mpc. We similarly fit £q 
in redshift intervals for the A cosmology (Fig. [l2]) using the 
60 ft -1 Mpc maximum as above. Again there is very little 
evidence of evolution. We note that there is some variation 
in the slope and amplitude of these power laws, but this 
appears to be mainly driven by the relatively low signal-to- 
noise in each redshift bin. Great care should be taken when 
trying to interpret these power law fit results, as amplitude 
and slope are correlated. 

An alternative method to derive a measurement of evo- 
lution is to constrain the power law slope and fit only for 



the scale length, so. This should be valid as we don't see any 
evidence for significant evolution in the slope of £q. We con- 
strain the slope to be that found over the full redshift range 
(Section^), 7 = 1.58 for the EdS cosmology and 7 = 1.56 
for the A cosmology. The results of this fitting process are 
seen in Table |l] and plotted in Fig. [l3|. In Table |l| we also 
list the reduced x 2 values for these fits. Limiting the fit 
to one parameter does not significantly alter the \ 2 values, 
demonstrating that the redshift averaged power law slope is 
a reasonable description of the data at all redshifts. Fig. [l3^ L 
shows that in the EdS cosmology clustering is constant as 
a function of redshift. The A cosmology result is shown in 
Fig. [l^b. In this case there appears to be a marginal increase 
by a factor of ~ 1.4 from z = 0.7 to 2.4. We compare our 
results in the EdS cosmology to previous QSO clustering re- 
sults from CS96 and LAC98, using their measurements of £ 
from 10 and 15 h~ x Mpc to obtain a value of so assuming 
a 7 = 1.58 power law (the best fit power law slope). Our 
results are in disagreement with those of LAC98 who find a 
~ 2a increase in clustering between z = 0.95 and z = 1.8. 
A possible cause of this is cosmic variance as LAC98 carry 
out their analysis in a single 24.6 deg 2 area of sky. However, 
given the large errors on the LAC98 data points, they only 
disagree with the 2QZ results at ~ 2a at 2 = 1.8. 

A non-evolving clustering distribution has strong impli- 
cations for models of structure and QSO formation. We first 
compare the 2QZ data to the simplest possible model, that 
of linear theory gravitational evolution in an Qo = 1 uni- 
verse. This model is applicable when QSOs either directly 
trace the mass distribution, or have a bias which is constant 
as a function of redshift. When fitting linear theory to the 
evolution in so for the EdS cosmology we find that the model 
is rejected by the 2QZ data at 99.8 per cent confidence. In 
the A cosmology the linear theory evolution rate is reduced. 
However in Fig. jjl^) we see that so increases with redshift, 
although the significance of the increase is marginal: a con- 
stant so as a function of redshift is not rejected by the data. 
When we try to fit linear evolution in this case it is rejected 
at > 99.9 per cent significance. If we require that the nor- 
malization of the mass clusterin g be fixed by eit her the local 
abundance of massive clusters ( Eke et al. 1996 ) or the 4-yr 
CO BE results ( Bennett et al. 1996] ) then the mass clustering 
scale length is forced to be less than so(z = 0) ~ 5 h -1 Mpc. 
In this case linear theory evolution is even more clearly re- 
jected by the 2QZ data. It therefore appears that QSO clus- 
tering cannot follow the linear evolution of the density field, 
and QSO bias must be a function of redshift. 

We should also make comparisons to galaxy clustering 
measurements. The typical scale length found in local galaxy 
surveys is so ~ 5 — 6 h -1 Mpc, only marginally higher than 
the 2QZ results for the EdS cosmology, and identical to the 
values found in the A cosmology. At z ~ 3 Adelberger et al. 
(1998) find a scale length of ro ~ 4 — 6 h~ 1 Mpc for Lyman- 
break galaxies, depending on the assumed cosmology. This 
again is very similar to the results derived from the 2QZ. 



4.2 Comparison to biased models of clustering 
evolution 

In the previous section we showed that for viable cosmologi- 
cal models, with evolution based on the gravitational growth 
of structure, QSOs do not simply trace the density fiuctu- 
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Figure 11. The two-point correlation function for 2QZ QSOs as a function of rcdshift for the EdS cosmology. Rcdshift increases, left 
to right and top to bottom. In each plot the solid line is the best fit power law on scales < 35 h^ 1 Mpc. The dotted line is the best fit 
to all the QSOs in the redshift range 0.3 < z < 2.9 and is shown to aid comparison between redshift intervals. The points without error 
bars at £(s) = 0.001 are where there are zero QSO pair counts in a bin. These points are fully taken into account in the fitting process. 



© 0000 RAS, MNRAS 000, 000-000 



12 S.M. Croom et al. 




100 



0.001 



10 

s (h _1 Mpc) 



Figure 12. The two-point correlation function for 2QZ QSOs as a function of rcdshift for the A cosmology. Rcdshift increases, left to 
right and top to bottom. In each plot the solid line is the best fit power law on scales < 60 h^ 1 Mpc. The dotted line is the best fit to 
all the QSOs in the redshift range 0.3 < z < 2.9 and is shown to aid comparison between redshift intervals. 
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Figure 13. The best fit values for so with fixed 7 as a function of redshift in the a) EdS and b) A cosmologies. The solid lines show 
the best fit li near theory evolution model in each case. The dotted lines show the best fit long-lived QSO biasing model (bias model 1) 
from Section 4.2.1. The dashed lines show the the models of Matarrese et al. (1997) (bias model 2) for different values of the minimum 
halo mass. The dot-dashed lines show the best fit empirical bias model (bias model 3). In a) the two dotted and dot-dashed lines are for 
COBE (top at z = 0) and cluster (bottom at z = 0) normalization. 



ations in the Universe. Therefore QSOs are related to the 
mass distribution via a redshift dependent bias. The form 
of this bias depends on the physical mechanisms of QSO 
formation. 

The question of QSO lifetimes can be linked to their 
clustering. If QSOs have lifetimes which are cosmologically 
long (~ a Hubble time), this would imply that QSOs are 
intrinsically rare. They could therefore be hig hly clustered. 
existin g in rare high peaks in the density field ( Efstathiou & 
Rees 19,88|), assuming that halo mass is the dominant factor 
in QSO formation. However, QSOs could form in a 'random' 
subset of less biased haloes, with the formation being driven 
by mechanisms other than mass, e.g. angular momentum. 

Alternatively, QSOs could have shorter lifetimes, of the 
order ~ 10 6 - 10 s years. There is mounting evidence for 
this, with the suggestion that most nearby galaxies appear 
to contain central supermassive black holes (e.g. Magorrian 
et al. 1998), so that most galaxies pass through a QSO/AGN 
phase. They could therefore be clustered in a similar man- 
ner to galaxies. However, even if all galaxies go through a 
QSO/AGN phase, it is possible that this phase picks out a 
particular time in the evolution of galaxies, e.g. epochs of 
major star formation or merging. QSO clustering evolution 
can potentially help us to distinguish between a number of 
possible QSO formation mechanisms. However, we should be 
wary of over interpreting models which do not include the 
uncertain physical mechanisms required for QSO formation. 

4-2.1 A long-lived QSO model 

The next simplest assumption, after assuming that bias does 
not evolve with redshift, is that QSOs are long lived (with 
ages of order a Hubble time). We assume that after forma- 
tion at some arbitrarily high redshift the subsequent evolu- 



tion of QSO clustering is governed purely by their motion 
within the gravitational pot ential pro duced by the density 
fluctuations in the Universe (Fry 1996). This then implies a 
bias which evolves as 



6(«) = l + (6(0)-l)6?(fio,A 0) 2). 



(9) 



We call this bias model 1. G(Q,o,\o, z) is the linear growth 
rate of density perturbations, which for an EdS cosmology 
is 1 + z. For the cosmological dependence of the growth 
rate we use the accurate fitting formula of Carroll, Press & 
Turner (1992), which is good to a few per cent (note that 
our G(Qo, Ao, z) is the full evolution term, and shouldn't 
be confused with the function of Carroll et al. which only 
contains the cosmological dependence). The biasing model 
of Eq. [] is also equivalent to QSOs forming in peaks of the 
density field above a constant threshold (CS96) . This model 
places certain limitations on the form of evolution. First, 
bias will tend to unity as time increases. Secondly, positive 
evolution (an increase in clustering) as redshift increases is 
not possible. This is because at most the bias only evolves as 
fast as G(f2o, Ao, z), cancelling out the growth in the density 
field. 

For comparison to the observed clustering we have nor- 
malized the mass e volution in two ways; using both local 
cluste r abundances jEko et al. 1996 ) and the 4-yr COBE re- 
sults (Bennett et al. 1996). We calculate so{z) for the mass 



assuming a CDM power spectrum with a shape parame- 
ter of r e ff = 0.25 (varying the shape parameter T c g only 
has an impact on the normalization when using the COBE 
data). In the EdS cosmology the so fits give 6(0) = 1.82^°'^ 
(1.62± ;^) for cluster (COBE) normalization. These corre- 
spond to erg ~ 1 for QSOs at z = (as for mass fluctuations 
is 0.52 and 0.65 for cluster and COBE normalization respec- 
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tively) , which is the same as the nominal erg ~ 1 value found 
for local galaxies. 

Schade, Boyle & Letawsky (2000) find that at low red- 
shift typical QSOs and AGN (where by typical we mean at 
or around the break in the luminosity function) have host 
galaxies that are remarkably similar to normal galaxies, ex- 
cept for a bias towards spheroid dominated galaxies. Ap- 
proximately 55 per cent of their sample had hosts which 
were best fit by a bulge-only model. Elliptical galaxi es are 



well kn own to b e more strongly clustered than spirals (Love- 
day et ; ,1. 1995 ) with a relative bias factor of 6 C , S ~ 1.9. Cor- 



recting for this morphological segregation gives an expected 
(78 = 1-2 — 1.3 for QSOs at low redshift, approximately in 
line with the above value. 

For the A cosmology a biased model of the form in Eq. 
H provides only a marginally adequate fit to the data (re- 
jected at 88%) with a best fit bias of 6(0) = 1.84±g;g|. This 
is because the model cannot reproduce the increase in clus- 
tering strength at high redshift visible in this cosmology. The 
hypothesis that QSOs have cosmologically long (~ Hubble 
time) lifetimes therefore appears unlikely in the A cosmol- 
ogy- 



4-2.2 More general models of biasing 

The above simple model of biasing can be extended in a 
number of ways. The most obvious is to remove the con- 
straint that objects formed at an arbitrarily high redshift, 
and allow objects to continue to form at lower redshift. The 
problem then becomes one of deciding how and when ob- 
jects do form. A natural method for deciding when dark 
matter haloes form is based on an application of the Press- 
Schechter (1974) formalism which describes the evolution of 
the number density of dark matter haloes. Working within 
this formalism Mo & White (1996) have obtained an approx- 
imation for the linear bias of dark matter haloes as a func- 
tion of mass. Matarrese et al. (1997) have used these ideas 
to provide biasing models in a CO BE normalized Qq = 1 
universe assuming a CDM power spectrum with a shape pa- 
rameter of r c ff — 0.25. These were extended to a number of 
different cosmological models by Moscardini et al. (1998). In 
particular we are interested in the transient model of Matar- 
rese et al., so called because the model does not require a 
normalization at z — 0. In this model, one assumes that all 
objects exceeding a given mass cut off can be observed at 
any given redshift. The bias (which we call model 2) then 
has the form 



b(z) = 1 - 1/Sc + [6(0) - (1 - l/tf c )]G(fio, Ao, zf 



(10) 



where 8 C is the critical linear overdensity for spherical col- 
lapse. For an EdS cosmology 5 C — 1.686 for all redshifts, 
however it only varies away from this value b y a few pe r 
cent for the other cosmologies considered here (Lilje 1992). 
Matarrese et al. find the values of 6(0) and by fitting to 
their Press-Schechter based models. These parameters de- 
pend on the minimum halo mass M m i n considered. We com- 
pare this model of biasing to QSO clustering is an EdS uni- 
verse in Fig. [Cm for minimum halo masses of M m i n = 10 11 , 
10 12 and lO x ^Vf0. In this cosmology the data are approxi- 
mately consistent with a minimum halo mass of 1O 12 M0 (al- 
though the model is still too steep) , while the normalization 
is too low (COBE normalization) for lower mass haloes, and 



the evolution is too steep for higher mass haloes. In the A 
cosmology, we compare the COBE normalized ACDM model 
of Moscardini et al. to our data. We note that this model is, 
in fact, for flo = 0.4, Ao = 0.6. However, given the model and 
data uncertainties these are adequate to make a general com- 
parison to the 2QZ clustering evolution in the A cosmology. 
In this case we find that the data are more consistent with 
(although slightly above) a model with M m i n ~ 1O 13 M0. 

Although these models appear to adequately describe 
the clustering evolution of QSOs, it is not at all clear what 
the physical justification for this is. The models of Matar- 
rese et al. assume that at each redshift QSOs inhabit the 
same mass haloes; this need not necessarily be the case. For 
example, Percival & Miller (1999) compare the evolution of 
bright QSOs, -25.4 > M B > -27.9, to the dark matter halo 
formation rate in a number of cosmologies. They find that 
for an EdS universe, with a CDM-type power spectrum of 
shape parameter T c g = 0.25 which is cluster abundance nor- 
malized, the evolution of bright QSOs is best fit by haloes 
of mass ~ lO la6 M . Our A cosmolo gy increases the mass 
to ~ lO n - 8 M . These masses are ~ 10 x smaller than those 
required to fit the 2QZ QSO clustering according to the 
models of Matarrese et al. This serves to demonstrate that 
we should be wary of over interpreting fits to models which 
do not contain a physical description of QSO formation. For 
example, it is possible that QSO clustering is a function of 
luminosity, a point which has not been discussed in this pa- 
per, but will be investigated in future work. 



4-2.3 An empirical biasing description 

Lastly we fit a purely empirical biasing model to the data. 
For this model we use a generalization of Eqs. ^ and [to] which 



b(z) = l+(b(0)-i)G(n ,\o,zf, 



(11) 



where 6(0) and /3 are left free to be determined by fitting to 
the data. We call this form of bias evolution model 3. The 
normalization of the mass density field is set by either cluster 
or COBE normalization as in model 1. The dot-dash lines in 
Fig.^ show the best fit empirical model for each of our as- 
sumed cosmologies. In the EdS case we find 6(0) = lA5t° il 
and (3 = 1.681°;** (6(0) = 1.28tg;J? and /3 = 1.89±g^) for 
cluster (COBE) normalization. As we might expect the A 
cosmology has a larger f3 with the best fit parameters be- 
ing 6(0) = 1.20±S:m and I 3 = 2 - 75 ±o.57- The relatively high 
normalization in this cosmology and the slow rate of mass 
clustering evolution means that a large value of j3 is required 
to fit the data. 



5 CONCLUSIONS 

The preliminary release dataset of the 2QZ contains 10681 
QSOs. It is already a factor of ~ 25 larger than previous 
QSO surveys to this depth (6j < 20.85). When completed 
the full sample will contain ~ 25000 QSOs. The current data 
set already allows us to measure the clustering of QSOs to 
un-precedented accuracy. In particular we find: 

1) QSO clustering integrated over the redshift inter- 
val 0.3 < z < 2.9 is well fit by a power law on scales 
~ 1 — 35 Mpc. In an Einstein-de Sitter universe the best 
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fit power law has s = 3.99±g:|| h' 1 Mpc and 7 = 1.581^. 
Introducing a cosmological constant increases the distances 
between QSOs, so that the scale length of clustering in- 
creases also. The power law then extends to ~ 60 hT 1 Mpc 
and is best fit by s = 5Wt :f h' 1 Mpc and 7 = 1.56i£j$. 
These results are remarkably similar to the clustering of nor- 
mal galaxies locally (z ~ 0.05). 

2) We compare the clustering of 2QZ QSOs to the 
ACDM model and find that the shapes of model and data 
are consistent. A comparison to a family of CDM models 
with different shape parameters, F c ff, finds that r c fr = 0.2 
to 0.4 provides an acceptable fit in the EdS cosmology. In 
the A cosmology only r e fl = 0.1 or 0.2 provide an acceptable 
fits due to the movement of structure to larger scales. This 
suggests a test for cosmological parameters using the linear 
break in the correlation function which will be possible using 
the completed 2QZ data set. 

3) We measure the clustering amplitude of QSOs as a 
function of redshift, parameterized by so assuming a fixed 
power law slope. In an Einstein-de Sitter universe we find 
that QSO clustering is constant in comoving coordinates 
over the entire redshift range we probe. In a A dominated 
universe we find that clustering appears to increase (al- 
though constant clustering is not excluded) with increasing 
redshift. For both EdS and A cosmologies a model in which 
QSOs follow the same evolution as linear theory gravita- 
tional clustering (or have a bias which is constant as a func- 
tion of redshift) is rejected at the > 99 per cent level. If the 
constant clustering is extrapolated to z ~ 3 it comfortably 
overlaps the clustering amplitu de found for Lyman-break 
galaxies (Adelberger et al. 199S). 

4) We compare simple redshift dependent bias mod- 
els to the measured clustering evolution. We first use a 
model in which QSOs are long lived (on cosmological time 
scales), so that their clustering simply evolves according to 
their motion in the gravitational potential. This is consistent 
with 2QZ clustering evolution in an EdS case, and predicts 
(Ts(z = 0) ~ 1 for QSOs, which is consistent with galaxy 
clustering. The long lived model is not able to reproduce 
the increase in clustering seen in the A cosmology and is 
marginally rejected at 88 per cent confidence. More com- 
plex models of QSO bias based on the Press-Schechter for- 
malism, have been developed by a number of authors. We 
use the models of Matarrese et al. (1997) and Moscardini et 
al (1998) to make comparisons to the evolution of the 2QZ 
data set. These models adequately describe the 2QZ clus- 
tering evolution when the minimum halo mass considered is 
M min ~ 1O 12 M (EdS) or M min ~ 10 13 M© (A). However, 
without a convincing model of QSO formation, the inter- 
pretation of the comparison to these models of clustering 
evolution is questionable. We lastly derive a fit to an empir- 
ical biasing model based on power law evolution of bias. 

The large volumes sampled by QSO surveys allow struc- 
ture to be investigated on the scales where growth is gov- 
erned by linear theory. Thus, meaningful measurements of 
large-scale structure, that are easily related to the underly- 
ing cosmology, can be made irrespective of the relative bias 
of QSOs. QSOs therefore play an crucial role in linking low- 
redshift/small-scale galaxy clustering measurements to the 
fluctuations in the density field at high redshift seen in the 
cosmic microwave background. The completed 2QZ survey, 
without the current varying observational coverage, will al- 



low detailed measurements of structure on a range of scales 
from ~ 1 to 1000 ft" 1 Mpc. 
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